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Abstract 



An existence result is proved for a nonlinear diffusion problem of phase-field type, 
consisting of a parabolic system of two partial differential equations, complemented 
5^ , by Neumann homogeneous boundary conditions and initial conditions. This system 

is meant to model two-species phase segregation on an atomic lattice under the 
presence of diffusion. A similar system has been recently introduced and analyzed 
in [3]. Both systems conform to the general theory developed in [5]: two parabolic 
PDEs, interpreted as balances of microforces and microenergy, are to be solved for 
the order parameter p and the chemical potential /x. In the system studied in this 
note, a phase-field equation in p fairly more general than in |3j is coupled with a 
highly nonlinear diffusion equation for p, in which the conductivity coefficient is 
allowed to depend nonlinearly on both variables. 
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Existence for a strongly coupled Cahn-Hilliard system 



1 Introduction 

In this paper, we prove an existence result for the following system in the unknown fields 
p and p: 



{l + 2g{p))dtfi + pg'{p)dtp-dw{K{p,p)Vfi) =0, (1.1) 

dtp-Ap + f\p)=pg'ip), (1.2) 

(^(/i, p)V/i) ■ z/|r = and 9^p|r = 0, (1.3) 

/i(0) = po and p(0) = po. (1.4) 



Each of the partial differential equations fll.ip - fll.2p is meant to hold in a three-dimensional 
bounded domain Q, endowed with a smooth boundary F, and in some time interval [0, T]. 
Such a system generalizes the phase-field model of Cahn-Hilliard type studied recently in 
[3]. Both models are of the type proposed in [5], and aim to describe phase segregation 
of two species (atoms and vacancies, say) on a lattice in presence of diffusion. The state 
variables are the order parameter p, interpreted as the volume density of one of the two 
species, and the chemical potential p. For physical reasons, p is required to be nonnega- 
tive, while the phase parameter p must, as such, obey < p < 1. Here are the features 
of [3] that have been generalized. 

Firstly, the nonlinearity / considered in [3] is a double-well potential defined in (0, 1), 
whose derivative /' diverges at the endpoints p = and p = 1: e.g., for / = /1-I-/2 with /2 
smooth, one can take /i(p) = c (p log(p) -|- (1 — p) log(l — p)), with c a positive constant. 
In this paper, we let fi be a maximal monotone graph from M to M. Consequently, 
equation (11.21) has to be read as a differential inclusion, in which the derivative of the 
convex part /i of / is replaced by the subdifferential /3 := dfi, i.e., 

dtp-Ap + ^ + f!,ip)=pg'ip) with ee/3(p); (1.5) 

moreover, since /i is not required to be smooth, its subdifferential may be multivalued; 
the selection of ^ in /3(p) is a further difficulty we face. 

Secondly, while in [3] g{p) = p, here g is any nonnegative- valued smooth function, 
defined (at least) in the domain where /i and its derivative (or rather, its subdifferential) 
live. 

Thirdly, and this is the most important novelty, conductivity k is not anymore a con- 
stant, but rather a positive- valued, continuous, bounded, and possibly nonlinear, function 
of p and p. For simplicity, we confine ourselves to study the existence of a solution under 
an assumption that guarantees uniform parabolicity, i.e., k > > 0. We point out that 
in a recent study [1] we let k depend only on p and possibly degenerate somewhere. 

Finally, relations (11. 4p specify the initial conditions for p and p, while (II. 3p are noth- 
ing but homogeneous boundary conditions of Neumann type, involving precisely those 
boundary operators that match the elliptic differential operators in (ll.ip - (ll.2p . 

Our paper is organized as follows. In the next section, we state our assumptions and 
our results. The existence of a solution is proved in Section [3l making use of a time-delay 
approximation and of a number of a priori estimates, that allow us to pass to the limit 
by compactness and monotonicity techniques. 
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2 Results 



In this section, we describe the mathematical problem under investigation, make our 
assumptions precise, and state our results. First of all, we assume 1] to be a bounded 
connected open set in with smooth boundary F (treating lower- dimensional cases would 
require only minor changes). Next, we fix a final time T G (0, +oo) and set: 

g:=fix(0,T), S:=Fx(0,T), (2.1) 
V ■= H\n), H := L'^{n), W := {v G H\n) : d^v = on F}. (2.2) 

We endow the spaces (12.21) with their standard norms, for which we use a self-explanatory 
notation like || ■ \\v', for powers of these spaces, norms are denoted by the same symbols. 
We remark that the embeddings W G V G H are compact, because Q is bounded and 
smooth. Moreover, for p G [l,+oo], we write || ■ ||p for the usual norm in L^^Q); as no 
confusion can arise, the symbol || ■ \\p is used for the norm in U'{Q) as well. 

First of all, we present the structural assumptions we make. We require that: 

K : {m,r) H- K{'m,r) is continuous from [0, +oo) x R to M, (2.3) 

the partial derivatives drK and d'^n exist and are continuous, (2.4) 

G (0,+oo), (2.5) 
< K,(m,r) < K,*, \drK{m,r)\ < k*, \d'^K{m,r)\ < k* for m > and r G M (2.6) 

/ = /i + /2, /i :M^ [0,+oo], /2:M^M, : M ^ [0, +oo), (2.7) 

/i is convex, proper, l.s.c. and /2 and g are functions, (2.8) 

/a, (?, and g' are Lipschitz continuous. (2.9) 

For convenience, we set: 

k' := drK, k" := d'^K, /3 := 9/i , and vr := /g; (2-10) 

pm pm pm 

K{m,r) := / K{s,r)ds, Ki{m,r) := / K'{s,r)ds, K2{m,r) := / K"{s,r)ds 
Jo Jo Jo 

for m > and r G M; (2.11) 

and we write D{f\) and -D(/3) for the effective domains of f\ and (3, respectively. Clearly, 
thanks to 

max{|i^(m, r)|, \Ki{m,r)\, \K2{m,r)\} < K*m for every m > and r G M. (2.12) 

We also note that the structural assumptions of [1] are fulfilled if k only depends on m, 
and that a strong singularity in equations (11. 2p for p is allowed. At variance with [1], 
equation (11. II) for /i is here uniformly parabolic, since g is nonnegative and k is bounded 
away from zero. 

Remark 2.1. Note that any convex, proper, l.s.c. function is bounded from below by 
an affine function (cf., e.g., [U Prop. 2.1, p. 51]), so that the assumption /i > looks 
reasonable, because one can suitably modify the smooth perturbation /2. Moreover, let 
us point out that the other positivity condition, (7 > 0, is just needed on the set -D(/3), 
while g can be extended outside of -D(/3) accordingly. 
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As to initial data, we require that: 

/io G 1/ n L'^in), PoEW, /io > and po e D{(3) a.e. in fi; (2.13) 
there exists E H such that ^ /^(Po) a.e. in f2. (2-14) 

Since /i is convex and /2 smooth, the above assumptions imply that /(po) £ 

As to the a priori regularity we require for any solution (/i, p, ^), we begin to observe 
that, for any given p, equation fll.Sp has the form of a standard phase-field equation. 
Therefore, it is natural to look for pairs (p, that satisfy 

p e W^'°°{0, T; H) n H\<d, T; V) H L°°(0, T; W), (2.15) 
eeL°°(0,T;i7). (2.16) 

Note that the Neumann boundary condition for p has been incorporated into fl2.15p 
(cf. fl2.2P o). Next, as to p, we require that 

l^e H\0,T;H)nL°°{0,T;V)nL'^{Q), p>0 a.e. in Q, (2.17) 
div(K(p,p)Vp) e L2(0,T;i7) and (k(p, p)Vp) ■ z/ = a.e. on S, (2.18) 

and note that we can expect that p G L^(0,T; W) (from the regularity theory of elliptic 
equations) only if the function k is smooth with respect to both variables. Nevertheless, 
(I2.17P and the regularity of the divergence are sufficient to write the Neumann boundary 
condition as done in fl2.18p . We also observe that 

peC\[0,T];C\n)) = C'(Q), (2.19) 

as a direct consequence of (12.151) and the compact embedding W C C°(fi) (see, e.g., [6l 
Sect. 8, Cor. 4]), whence g'{p) G C°(Q). Thus, under all of the above requirements, we 
can write the system of equations and the initial condition in the following strong form 

{l + 2g{p))dtiJ + i2g'{p)dtp-dw{K{i2,p)Vi2) =0 a.e. in Q, (2.20) 
dtp- Ap + ^ + n{p) = pg'{p) and ^ G /3(p) a.e. in Q, (2.21) 
p(0) = Po and p(0) = po a.e. in Q. (2.22) 

Here is our existence result. 

Theorem 2.2. Assume that both flOj) - fl23TD and i^A^-(^J^ hold. Then, there exists 
at least a triplet (p,p, ^) satisfying (12. 15p - (12.181) and solving problem (I2.20l) - fl2.22p . 

This is the only result we prove in the present paper. We note, however, that the 
uniqueness result obtained in [1] still holds here, provided that n is taken constant and 
Po smoother. For the reader's convenience, we summarize the results of [1]. 

Theorem 2.3. Assume that both flO)- (E:TT]) and (I^TT^ - (Oip hold, and, moreover, 
that Po G ly and k = 1. Then, there is a unique triplet (p,p, ^) satisfying (I2.15p - (l2.18p 
and solving problem (12. 20p - (12.220 and its component p enjoys the following regularity 
property: 

/i G H^^'P(0,T;i7) nLP(0,T;iy) for every p E [1, +oo). (2.23) 
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Throughout the paper, we make use of some well-known embeddings of Sobolev type, 
namely, V C L^{Q) for p e [1, 6], together with the related Sobolev inequality 

IK'IIp ^ C'll'^llv for every v & V and 1 < p < Q, (2.24) 

and W^'P{Q) C C^{Tl) for p>3, together with 

IK'lloo < C'p IK' II VKi>p(n) for every v G W^'^^Q) and p > 3. (2.25) 

In fl2.24p . C depends only on Q, while Cp in f l2.25p depends also on p. In particular, the 
continuous embedding W C W^'^{fl) C C°(fi) holds. Some of the previous embeddings 
are in fact compact. This is the case for V C L"^(r2) and W C C^{Q). We also account 
for the corresponding inequality 

IK'IU < £^||Vi;||/f + Ce|K;||// for every e and £> (2.26) 

where depends on Q and e, only. Furthermore, we repeatedly make use of the notation 

Qt:=nx (0, t) for t E [0, T], (2.27) 

of the well-known Holder inequality, and of the elementary Young inequality 

ab <ea^ + ^ iP' for every a, 6 > and e > 0. (2.28) 

Finally, again throughout the paper, we use a small-case italic c for different constants, 
that may only depend on f2, the final time T, the shape of the nonlinearities / and and 
the properties of the data involved in the statements at hand; a notation like Cg signals 
a constant that depends also on the parameter e. The reader should keep in mind that 
the meaning of c and might change from line to line and even in the same chain of 
inequalities, whereas those constants we need to refer to are always denoted by capital 
letters, just like C in (EMj). 



3 Existence 

In this section, we prove Theorem [221 which ensures the existence of a solution. Although 
our proof follows the argument in |3] and |1] closely, we present the whole argument, and 
sometimes give some details, since the changes with respect to the quoted papers are 
spread over the whole exposition. Our starting point is an approximating problem, which 
is still based on introducing a time delay in the right-hand side of fl2.2ip . Precisely, we 
define the translation operator 7j. : L^(0,T; H) — > L^(0,T; H) depending on a time step 
r > by setting, for v G L^(0, T; H) and for a. a. t G (0, T), 

i7rv){t) := v{t - r) if t > r and {7rv){t) := /iq if t < r (3.1) 

(the same notation 7t-v will be used also for a function v that is defined in some subinterval 
[0, T'] of [0, T]). At bottom, what we do is to replace /i by T,-// in fl2.2ip . However, since it 
is not obvious that we can keep /i positive, we extend k to a function k : R — t- M satisfying 
similar properties. Moreover, we assume that the analogue of (12.60 holds for R and its 
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derivatives, with the same constants and k* (we replace and k* by 2^^ and k*/2 in 
the original (12.61) if necessary). So, the approximating problem consists of the equations 

{l + 2g{p^))dtHr + lJ'Tg'{pT)dtpr-div{R{fi-r,pT)'^fJ'T) =0 a.e. in Q, (3.2) 
dtPr - Apr + ^r + t^{Pt) = {^tPt) q' {Pt) and € /3(pr) a.e. in Q, (3.3) 

complemented by the initial and boundary conditions 

/i^(0) = /io , Pr(0) = Po , dyPr\T. = 0, (^(/ir, PtW Pt " = 0. (3.4) 

For convenience, we allow r to take just discrete values, namely, r = T /N , where is 
any positive integer. 

Lemma 3.1. The approximating problem has a solution [fir, pr, C,t) satisfying the analogue 

of (^M-<^M- 

Proof. We confine ourselves to give a sketch. As in [3], we inductively solve problems 
on the time intervals /„ = [0,t„] := [0,r;,r], n = 1,. . . ,N, by constructing the solution 
directly on the whole of at each step. Namely, given Pn-i, which is defined in f2 x /„_i, 
we note that T^/i^.i is well defined and known in f2 x and solve the boundary value 
problem for p„ given by the phase field equations 

dtpn- ^Pn + ^n + 'n-iPn) = {7rPn)g'{Pn) and ^„ G /3(pn) in X J„ (3.5) 

complemented by the boundary and initial conditions just mentioned for p^-. Such a 
problem is quite standard and has a unique solution p„ in a proper functional framework. 
Now, we observe that the function 

k : {x,t,m) R{m, pn{x,t)), {x,t) G x /„, m G M 

is a Caratheodory function satisfying k,^ < k < k* on its domain, so that the equation 

(l + 2g{pn)) dtpn + Png'iPn) dtPn - div (k(p„, p„) Vpn) =0 m Q X In (3.6) 

in the unknown function p„ is uniformly parabolic (let also recall that g is nonnegative). 
Thus, we can solve the problem obtained by complementing (13. 6p with the boundary and 
initial conditions prescribed for pr- Therefore, the problem to be solved has a unique 
solution in a proper space, provided that the coefficient g'{pn)dtPn is not too irregular. 
So, we should prove that, step by step, we get the right regularity for p„ and This 
could be done by induction, as in [5], with some modifications due to our more general 
framework. We omit this detail and just observe that the needed a priori estimates are 
close (and even simpler, since r is fixed here) to the ones we perform later on in order 
to let r go to zero. The final point is p„ > 0. We give the proof in detail. We multiply 
equation (13.61) by — p^ := — (— p^)"*", the negative part of p^, and integrate over Qt with 
any t ^ I^.. We observe that 

[{l + 2g{pn{t)))dtPn + Png{Pn)dtPn\ {- P') = ^ {{1 + 2 g (pn)) \ P^ \^) ■ 

Hence, by using po > and owing to the boundary condition, we have 

^ / (1 + 2(7(p„(t))) |p;(t)|' + / k{pn,Pn)\'^Pn\^ = for every t G /„. 
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As both g and R are nonnegative, this imphes /i^ = 0, that is, /in > a.e. in x J„. Once 
all this is checked, the finite sequence {fin, Pn,^n), n = 1,. . . ,N, is constructed and it is 
clear that a solution to the approximating problem we are looking for is simply obained 
by taking n = N. □ 



Thus, we fix a solution for each r. We note that, a posteriori, we can 

replace k by k in (13. 2p . since the component fir of our solution is nonnegative. Our aim 
is to let r go to zero, so as to obtain a solution as stated in Theorem 12.21 Our proof 
uses compactness arguments and thus relies on a number of uniform (with respect to r) 
a priori estimates. In order to make the formulas to come more readable, we write fi and 
p rather than fij. and Pr in the calculations. 

Remark 3.2. Sometimes, when deriving our a priori estimates, we proceed formally. 
However, our procedures can be made rigorous. For instance, one can get more regularity 
for the approximating problem by regularizing k and the initial data, if necessary. More- 
over, local regularization is often sufficient. Consider, e.g., equation (13. 3p and rewrite it 
in the form 

-Ap + p + /3{p) 3h:=p-dtp- 7i{p) + {7r^Ag'{p). (3.7) 
Now, (here t is just a parameter) the elliptic equation: 

— An + u + w = h and w G /3(n), 

complemented by homogeneous Neumann boundary condition, yields a well-posed prob- 
lem; as is well known, its solution (m, w) is the limit of the much smoother pair (m^, /^^(Me)), 
where fig is a regularization of {3 of Yosida type (see, e.g., [21 p. 28]; see also the proof 
of Lemma 3.1 of [3] for a further regularization) and is the solution of the analogous 
boundary value problem for 

- AUe + Ue + /3e{Ue) = h. 

On the other hand, we have {u,w) = (p, by (13.71) . Therefore, it is essentially correct 
to regard /3 as if it were a smooth function in the original equation (13. 3p . and treat such 
equation like a more regular one (e.g., by differentiating it or taking irregular functions 
as test functions). 

First a priori estimate. We test (13. 2p by Pr and observe that 

[(1 + 2g{p)) dtp + p g\p) dtp] p = ]^dt[{l + 2g{p))p''] . 

Thus, by integrating over (0,t), where t G [0,T] is arbitrary, we obtain 

/ {l + 2g{p{t)))\p{t)\^+ [ ^{p,p)\Vp\'= [ {l + 2g{po))pl. 
Jn JQt Jn 

Hence, we recall that g >0 and ft > > 0, and conclude that 



Pr||L°°(0,r;/f) + ||Pr||L2(0,T;y) < C. 



(3.8) 
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Consequence. The Sobolev inequality ( I2.24p . estimate (I3.8p . and (I2.12p . imply that 

||Atr|U2(o,T;L6(n)) + || ^ (/ir , Pr ) || (0,T;H)nL2(o,T;L6(Q)) < C with ijj = K, Ki, K2. (3.9) 

Another implication of (13. Sp . along with (13. ip and (I2.13p . is 

\\7rlJ.r\\L°°iO,T;H) + \\'^tIJ't\\l^(0,T;V) < C. (3.10) 

Second a priori estimate. We add Pr to both sides of (13.31) and test by dtPr- We 
obtain: 

/ M + l\\p{t)ry+ I Mp{t)) 
JQt ^ Jn 

= \l |Vpor+ / f{po) + \ I {p\t) + 2h{p{t)))+ I g\p){7rli)dtp 

<C + C [ |p(t)P + i / M + c\\7rl^\\lo.^0,T;H), 

Jn JQt 

for every t G [0,T]. Thanks to the chain rule and the Young inequality (12. 28 p . we have: 



1 



c I \pit)\'<c \po\' + - / M + c / \\pis)\\lds. 



n Jn ^ JQt Jo 

Then, as /i is nonnegative, by accounting for (13.80 . with the help of the Gronwall lemma 
we infer that 

' \dtp\' + ymi + [ fMt))<c. 



Jn 

Thus, we conclude that 

||Pr||//i(0,r;H)nL-(0,r;y) < C and \\fiPr)\\L^(O^T;L^n))<C. (3.11) 

Third a priori estimate. We proceed formally (see Remark 13.20 . We rewrite (13.30 as 

-Ap + (3{p) = h:= -dtp - 7T{p) + {7,fi)g'{p), (3.12) 

and multiply this relation by either —Ap or /3(p). By doing that, we derive an estimate 
for both Ap and /3(p) and we can use the regularity theory for elliptic equations. We 
conclude that 

\\Pt\\lHO,T;W) < C and Ur\\L2(0,T;H) < c (3.13) 

Fourth a priori estimate. Our aim is to improve estimates (13. lip and (13.131) . To 

do that, we proceed formally, at least at the beginning (our procedure could be made 
completely rigorous, as sketched in Remark l3.2p . We start from an estimate coming from 
the theory of maximal monotone operators [2], namely, 

||5iM(0)||H< 11^(0) + Apo||h+ min llr^lln, (3.14) 
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for the unique solution {u,u) to the equations (cf. (13. 3p ) 

dtu — Au + u = ip := g'{p)7rfi — 7r(p) and u G f3{u), 

complemented by the same initial and boundary conditions prescribed for p. Note that 
in (I3.14P P is understood as the induced maximal monotone operator from H to H. By 
observing that (m, u) = (p,^), applying ( I3.14p . and combining with our assumptions on po 
(see (I2.14p . in particular), we obtain: 

\\dtPrm\H < c{\\fio\\H + llpolk + 1 + Uo\\h)= c. (3.15) 

We use (I3.15P in the calculation we are about to start: once again, we proceed formally, 
and write ^ = /3(p) as if /3 were a smooth function. We differentiate (13. 3 p with respect to 
time and test the equation obtained for dtp. We find: 



Qt 



I [ \dtp{t)\'+ [ mp\'+ [ p'ipMpi' 

^ Jn JQt JQt 

\{dtpm\'- [ n'{pMp\'+ [ g"{p){7^f^Mp\' 

g'{p)dt{7rP')dtp 

't 

<\ I mpm\' + c[ {l+'JrP)\dtp\'+ [ g\p)^t{7r^i)^tp. (3.16) 
^ Jn JQt JQt 

We treat each term on the right-hand side, separately. The first one is estimated by (13.150 . 
In order to deal with the second one, we account for the Holder inequality, (I3.10p . the 
compact embedding V C L^{VL) (see (12.260 ). and (13.110 . We obtain: 

/ {l+7rPi)\dtP? < f \\l + {7rP^){s)\\H\\dtp{s)\\lds 
JQt JO 

<c f \\dtp{s)\\lds<s f \Vdtp\^ + cJ \dtp\^ 
JO JQt JQt 

<e I |V9tp|2 + Ce , (3.17) 

JQt 

for every e > 0. Let us come to the last term of (13.161) . Firstly, on recalling that dt^^-fi = 
in (0, r) by the definition of 7^-, we compute dtfi from (13.21) . Then, we integrate by parts 
and have repeated recourse to Holder, Sobolev, and Young inequalities. We deduce that 



g\p)dt{7'rij)dtp= / / dtfi{s) dtg{p{s + t)) ds 
Jo Jn 



Qt Jo Jn 

Jn l + 2g{p{s)) [^^''('^^^^^)' ^^^))^'^^^)) ~ f^i^)9'ipis))dtp{s)] dtg{p{s + r)) ds 



t-T 



Jn 

rt-T 



.(M3),p(.))VM^)-vf^^^^rf. 

1 + 2^(p(s)) 



'0 Jn 



"''^!l'ff';;"" M^)a.P(..)a.P(.^ + r)d.. (3^8) 

1 + 2^(p(s)) 
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We treat the last two integrals separately, by using our structural assumptions. 
As to the former, we have: 

"t-T 



ri .(M.),p(s))vM3)-vf^M^c/. 

Jo Jn 1 + 2^(p(s)) 

's + r))dt 
1 + 2g{p{s)) 



Jn 



<c f [ \Vfi{s)\\Vdtp{s + T)\ds 
Jo Jn 

+ c / 7 \Wfiis)\\\/pis)\\dtpis + T)\ds 
Jo Jn 



'0 Jn 

+ c[ [ \Vpis)\\Vpis + T)\\dtpis + T)\ds. (3.19) 
Jo Jn 



Moreover, thanks to (13.81) . we infer: 



[ \Vfi{s)\\Vdtp{s + T)\ds<e [ \Vdtp\^ + cJ |V/ip 

Jn JQt JQt 



< e / \Vdtp\^ + c, , (3.20) 

JQt 

for every e G (0, 1). On the other hand, we also have: 

* 7 \Vpis)\\Vpis)\\dtp{s + T)\ds 
Jn 

< [ ^\\Wfiis)h\\\/pis)\\Mpis + r)\Uds 
Jo 

<e [\\d,p{s)ryds + cJ' ^||VM^)ll?.l|Vp(.)||^ci. 



<e[ \Vdtp\' + cf \dtp\' + cJ ^\\Vpi{s)\\W^p{s)\\lds 
jQt JQt Jo 

<el \Vdtp\' + c + Ce [ \\Vfi{s)\\U^p{s)\\lds (3.21) 
Jqi Jo 

(in the last inequality, (13. lip has been used). Now, we improve the estimate just obtained 
by owing to the regularity theory for linear elliptic equations, as well as to estimates (13. 8p 
and (13. lip . For each fixed s G (0,T), we have 

\\'^pis)rv<c{\\p{s)ry+\\Ap{s)\\i) 

<c + c\\-dtpis) - nipis))+g'ipis))7Ms)fH < \\dtp{s)\\j, + c. 
Therefore, the above estimate becomes 

Ct-T 



Jn 



\Vp{s)\\Vp{s)\\dtp{s + T)\ds 

<e [ \Vdtp\^ + c, [ \\Vfi{s)\\l\\dtp{s)\\lds + Ce. (3.22) 
Jot Jo 
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Analogously, one shows that 



t-T 

' \\/flis)\\Vpis + T)\\dtpis + T)\ds 

Jn 

t 



<e [ \Vdtp\^ + cJ \\V{7rP){s)\\l\\dtp{s)\\lds + c,. (3.23) 

JQt Jo 

Thus, by collecting ( K20\f and fl3:22|) -f l3:23|) . we deduce that f l339|) yields: 

<e [ \Vd,p\' + cJ (||VMs)||?,+ ||V(T./i)(s)||^)||9,p(s)||^d. + c,, (3.24) 
for every e > 0. 

We now take up the last integral in (13.181) . By using the compactness inequality (I2.26P 



and (13.81) . we have: 

^ g'{p{sy)g\p{s^r)) 
1 + 2g{p{s)) 



p{s)dtp{s)dtp{s + r) ds 



'0 Jn 

<C [ ^Ms)\\4dtp{s + T)\\4dtp{s)hds 







t-T ft 

2 J„ , „ / II ,./„M|2| 



<e / WdAs + rWyds + c, / Ms)rMp{s)rHds 
Jo Jo 

<e [ \Vdtpf + c + Cs [ Ms)\\l\\dtp{s)\\lds. (3.25) 

jQt Jo 

Therefore, due to (13:2^1) and (13:251) . fl3A8|) becomes: 

/ g\p)dt{7rP)dtP<2e f \Vdtp\'' + c, f \\V p{s)\\l\\dtp{s)\\l ds 

jQt JQt Jo 

+ Ce [ \\pis)\\l\\dtp{s)\\lds + Cs. (3.26) 
Jo 

At this point, we combine (13.151) . (I3.17p . and (I3.26P with (13.161) . and we choose e small 
enough. Since the last integral on the left-hand side of (I3.16P is nonnegative, because /i 
is convex, we obtain: 



/ \d,p{t)\'+ [ \Vdtp\'<c f<P{smp{s)\\lds + c, 
Jn Jot Jo 



where 0(.) := || V/i(s) ||^ + || V(T./i)(s) ||^ + ||/.(s) ||t 
As </) G L^(0,T) by (I3.8I) - (I3.10I) . we can apply the Gronwall lemma and conclude that 

\\dtPT\\L°°(0,T;H)nL'2{0,T;V) < C. (3.27) 

Consequence. Note that — Ap^ + = —dtPr + 9'{Pt)'^tPt is bounded in L°°(0, T; iJ), 
due to (13. 8p and (13.271) . Therefore, by a standard argument (multiply formally by ^r)^ 
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we deduce that both — Ap,- and are bounded in the same space, whence by elhptic 
regularity 

||Pr||L°°{o,T;iy) < c and ||Ct||l°°(o,T;H) < c; (3.28) 

moreover, 

||Pr||L-(Q) + ||^/'(Pr))||L-(Q) < c With tjj = g , g' , IT , (3.29) 

due to the continuous embedding W C L°°{Q) and the continuity of such ip^s. 

Fifth a priori estimate. To prove an L°° estimate rather than just a boundedness 
property, we borrow the argument in [3]. We observe that the approximating solution 
satisfies: 

^d,[{l + 2g{p))\{f,-k)+\'] = [{l + 2g{p))^,f, + {^,-k)g'{p)^tp]{^^-k)+, (3.30) 

for every k E M.. Hence, by assuming that k > yUg := ||uo||oo and by testing (13.21) with 
[fi — k)~^ , we obtain: 

\ ! {i + 2g{pmmt)-kY?^ j K{p,p)\v{fi-ky\' = -k [ dMif^-k^; 

^ Jn JQt JQt 

with this, on recalling that g >0 and k, > k^, we deduce the inequality 

2, \{p{t)-kY\^ + K. f \Vip-k)+\' = -k [ dtgip)ip-k)+. 
In [3], for e = 1, we have that g{r) = r and k = 1; the corresponding inequality is: 

^ [ m)-kr\'+ f \v{f,-ky\' = -k [ d,p{p-k)\ (3.31) 

^ Jn JQt JQt 

Therefore, the argument used in that paper can be repeated here essentially without 
changes. As a matter of fact, the analogue of (12.211) is never used in [5], the whole proof 
being based just on (I3.3ip . the regularity dtp G L°°(0,T; H) fl L^(0,T; V), and an upper 
bound, say Co, for the corresponding norm; moreover, the upper bound for p, that is 
constructed explicitely, depends only on fl, T, Pq, and Cq. In the present case, we have 
to use the same regularity for dtg{p) and estimate (13. 29 p . In conclusion, we obtain: 

||/-iT||L°°{Q) < c. (3.32) 

Sixth a priori estimate. We proceed formally, as done for the third a priori estimate, 
by writing ^ = /3(p) as if /3 were a smooth function (see Remark 13. 2p . We test by {^{t)Y 
(13. 3p . written in the form (I3.12p . at (almost) any fixed time t G (0,T). We obtain: 

Jn Jn Jn 

As the first term on the left-hand side is nonnegative, by the Holder inequality we de- 
duce that 



\\m\\i<\\Ht)\\e\\m)%/^ = \\Ht)h\mr 



6' 
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whence he have immediately that ||^(t)||6 ^ II ^(^) lie- We infer that 

l|Ap(t)||L6(n) < c||/i(t)||L6(n) and \\p{t)\\w^.6(^n) < c||/i(t) ||i6(n), 

first by comparison in fl3.12p and then by the standard regularity theory of linear elliptic 
equations. As W^'^{Q) is continuously embedded in C^{Q) (see fl2.25p ). and as the above 
inequalities hold for a. a. t G (0,T), we deduce that 

II Vp||L2(0,T;L°°(n)) < c||/i||l2(o,T;L6(q)). 

Now, we observe that h is bounded in L'^{0,T; L%Q)), thanks to flXTTj) . (IX^ . (K^ . 
and the Sobolev inequality. Therefore, we conclude that 

||Vpr|U2{0,T;L°°(n)) < c (3.33) 

A byproduct of our proof is that \\pT\\L^(o,T;W'^'(^{n)) < c. 



Seventh a priori estimate. On recalling f l2.10p - fl2.lip . the following preparatory iden- 
tities hold for the approximating solution: 

VK{fi, p) = k{p, p)V/i + Ki{p, p)Vp, (3.34) 

dtK{p, p) = k{p, p)dtp + i^i(/i, p)dtp, (3.35) 

dtKi{ii, p) = k\p, p)dtp + K2{p, p)dtp. (3.36) 

Moreover, we notice that 

\\^{iJ',p)\\l^{q) < c with ip = k,k',K,Ki,ot K2, (3.37) 

due to our structural assumptions 02.61) and to fl3.32p . f|2.12p . Now, we formally test f l3.2p 
by dtK{p,p) and get: 



[l + 2g{p))dtpdtK{fi,p)+ / pdtg{p)dtK{p,p) 

jQt 

+ [ K{p,p)Vfi-VdtK{p,p) = 0. (3.38) 
jQt 

With the help of fl3.34l) - fl3.35l) . we rewrite the first two integral as follows: 

:i + 2g{p))dtpm{p,p) 

:i + 2gip))n{fi,p)\dtp\'+ I {l + 2g{p))K,{p,p)dtpdtp 



Qt 



pdtg{p)dtK{fi,p) = / pdtg{p)K{p,p)dtp+ / fidtg{p) Ki{p, p) dtp 



Qt 
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In the third integral of f l3.38p . we also integrate by parts and use fl3.36p . We get: 

JQt JQt 

\VK{p{t),pm'- [ KMt),p{t))Vp{t)-VK{p{t),p{t))-c 
Jn 

VK{p,p)-dt{K^{p,p)Vp) 



'Qt 
1 



I I \VK{p{t),pm'- / KMt),p{t))Vp{t)-VK{p{t),p{t))-c 



+ / VK{p, p) ■ Vp {k\p, p)dtp + i^2(/i, p)dtp) + / p)VK{p, p) ■ Vdtp. 

JQt JQt 

If we insert these identities in (13.381) . on keeping just the positive terms on the left-hand 
side, on recalling that g > Q and k > k*, and on and accounting for estimates fl3.32p 
and fl3.37l) . then we deduce that 

/ M + \ I \VK{p{t),p{t))\^ 
JQt ^ Jn 

<c[ \dtii\\dtp\+c ! \dtp\' + c [ \Vp{t)\\VK{p{t),p{t))\ 
JQt JQt Jn 

+ c[ \VKip,p)\\Vp\{\dtp\ + \dtp\)+c [ \VKip,p)\\Vdtp\+c. (3.39) 
JQt JQt 

As the first three terms on the right-hand side can be trivially dealt with by accounting 
for (13.81) . (13.111) . and the elementary Young inequahty, we concentrate on the last two 
integrals. For every e > 0, we deduce that 

\VKip,p)\ |Vp| {\dtfi\ + \dtp\) + [ \VKip,p)\ \Vdtp\ 



Qt 



<e / \dt^i\' + ce / \S/p\'\WK{ii,p)[ 

1 



\VK{p,p)\^+^- I \Vp\'\dtp\'+^- I \Vdtp\ 
On the other hand, we have: 



Qt 2 JQ^ Z, JQ^ 



f \Vp\'\VK{p,p)\'+ [ |VA'(/i,p)p+ / \Vp\'\d,p\'+ [ mp\' 

JQt JQt JQt JQt 

< f<P{s)\\VK{p{s),p{s))\\lds+ f\\Vp{s)\\l\\dtp{s)\\lds+ [ mp\' 

Jo Jo JQt 

< [ (j){s) \\VK{p{s),p{s))\\lds+ ||Vp||i2(o,r;L-(n))ll^tPlli-(0,T;i?) + l|V5tp||i2(o,r;H). 

Jo 

where 0(s) := ||Vp(s)||^ -|- 1. As G L^{0,T) thanks to (13.330 . and as the last norms in 
the above inequality are bounded by (I3.27P and (I3.33p . we can choose e small enough and 
apply the Gronwall lemma. We conclude that 

\\dtlJ'T\\L^O,T;H) + \\K{Pr,pT)\\L--(0,T;V) < C. (3.40) 
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Consequence. By combining ( I3.40p . f l3.34p . and k > k*, we derive that 

||V/ir||L°°{o,T;H) < c, whence ||yU^||L<-{o,r;V) < c. (3.41) 
Furthermore, by comparison in f l3.2p . we also deduce that 

||div(K(/i^,p^)V/i^)||L2(o,T;H) < c. (3.42) 

Limit and conclusion. By the above estimates, there are a triplet (//, p, ^), with p > 
a.e. in Q, and two functions k and ( such that 

/i^^/i weakly star in i/i(0,T;i7) nL2(0,T;1/) nL°°(g), (3.43) 

pr^P weakly star in L°°(0, T; W^), (3.44) 

dtpr dtp weakly star in L°°{0,T; H) n L'^{0,T;V), (3.45) 

^ ^ weakly star in L°°(0, T; //), (3.46) 

K{fir, pr) k weakly star in L°°(0, T; V), (3.47) 

div(K(/iT-, p^)VPt-) -> C weakly in L^(0, T; if), (3 



at least for a susequence r = rj\0. By the weak convergence of time derivatives, the 
Cauchy conditions ( I2.22p hold for the limit pair (p,p). By f l3.43p -( l3.45p and the compact 
embeddings W C C°(0) and 1^ C if, we can apply well-known strong compactness results 
(see, e.g., [6l Sect. 8, Cor. 4]) and, possibly taking another subsequence, we have that 

/i^ — J- p strongly in L^(0,T; H) and a.e. in Q (3.49) 
pr^p strongly in C^(Q). (3.50) 

The weak convergence fl3.46p . together with (13.501) . implies that ^ G /3(p) a.e. in Q (see, 
e.g., [21 Prop. 2.5, p. 27]), due to the maximal monotonicity of the operator induced by (3 
on Li^iQ). Furthermore, the convergence stated in fl3. 490 - 03. 50p entails that 0(pr) — 0(p) 
uniformly in Q for (p = g, g', ir and ipifJ^r, Pt) — ^ i'ipi p) a.e. in Q for ip = k,K, whence, in 
particular, k = K{fi,p). As all the above functions tp{fir,PT) are uniformly bounded, we 
deduce that the convergence is in fact strong in U'{Q) for every p < +oo and weak star 
in L°°{Q) in each case. This shows that the limits of the products 

{1 + 2g{pr))dtfir, firg'{PT)dtPr, k(p^, p^) Vp^, aud firg'iPr) 

that appear in equations fl3.2p - fl3.3l) can be identified as the products of the corresponding 
limits. In particular, by using also 03.480 . we derive that div(/t(p, p) Vp) equals ( and 
belongs to L'^{Q). All this implies both O2.20p for the limit (p, p) and the convergence of 
the normal trace ^(pr, Pt) Vpr ■ i^- Thus, the expected Neumann condition also holds in 
the limit, and the proof is complete. 
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